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ABSTRACT 

35974 

After adjusting the estimate for the mass of a zero-absolute -visual- 
magnitude meteor, formulas for meteoroid mass are derived from E. J. Opik's 
( 1958) physical theory for dustball meteors. Masses are calculated for G. S. 
Hawkins and R. B. Southworth's ( 1958) random sample of 285 sporadic photo- 
graphic meteors. Weighting functions are derived for the statistical minimization 
of the effects of several sources of bias. Means, standard deviations, and ordi- 
nary, multiple, and partial correlations are shown for a 5-variate statistical 
analysis. Weighted cumulative distributions for parameters are plotted for two 
gradations with respect to mass. Data-point distributions are plotted which 
establish cumulative flux as functions of mass, momentum, kinetic energy, etc. , 
which show the effects of the earth’s motion on the distributions and correlation 
of parameters. 

The results support revisions in the author's (Jan. 1965) interpretation 
of Explorer XVI data and predicted puncture rates for the Pegasus meteoroid 
measurement satellites. The apparent discrepancy between the logarithmically 
linear relations between cumulative flux and mass, as inferred from photographic 
meteors by G. S. Hawkins and E. K. L. Upton ( 1958) and from zodiacal light and 
solar F-corona data by D. B. Beard ( 1963) is shown to be accountable by E. J. 
Opik's ( 1951) probability that a particular meteoroid will encounter the 
one revolution of the particle. 
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P 
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k , k 
Pi P2 

G o (v) 

G t <M p* 

n 
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Definition 

Meteor absolute visual magnitude. 

Meteor absolute photographic magnitude. 

Meteoroid mass in grams in space. 

Meteoroid velocity before deceleration in the atmosphere. 
Meteor total luminous efficiency. 

Effective visible length of a meteor trajectory. 

Effective time of meteor visibility. 

Luminous efficiency of meteor "initial" and "impact radiation. 

Luminous efficiency of meteor "temperature radiation." 

Value of m which would be indicated if (3 were replaced by /3 q . 

Meteor luminous efficiency for a diluted coma. 

Luminous efficiency for a compact cloud of meteor vapors. 

Degree of dilution of a meteor coma. 

Constants in the reduction formula for visual meteors. 

Constants in the reduction formula for photographic meteors. 

Functional representation of -log ( 3 ^. 

v) Functional representation of -log ( i + 13,/ (3 ) . 

t o 

Partial derivative of log (3 with respect to log v. 

Meteor height in kilometers above sea level at the point of 
maximum brilliance. 


t? 
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f b 
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6 
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Definition 

Meteor height in kilometers above sea level at the point of 
appearance. 

Meteor height in kilometers above sea level at the point of 
disappearance. 

Zenith angle to meteor radiant in radians. 

Elongation of the true radiant from the apex of the earth's 
way in degrees. 

Whipple's cosmic weight for meteors. 


Opik's meteoroid- earth encounter probability. 

Meteor heliocentric orbital eccentricity, or base of natural 
logarithms, depending on context. 

Celestial latitude of the corrected radiant in degrees. 

Weighting factor due to bias in height and velocity. 


• • 

Weighting factor for Opik's earth-encounter probability P. 

Right ascension of the corrected radiant in degrees. 

Declination of the corrected radiant in degrees. 

Limited reciprocal of the apparent fraction of the circle of 
celestial latitude through the corrected radiant. 

Difference between the celestial latitudes of the horizon points 
on the circle of celestial latitude /3 through' the corrected radiant, 
in radians. 

Celestial latitude of the zenith. 

Right ascension of the zenith in radians. 


viii 



Symbol 


6 

z 

t 

mu 


t 


du 


t 

s 


k 

z 


f 

o 




o 


Xi 


t • •• 



DEFINITION OF SYMBOLS (Cont'd) 

Definition 

Declination of the zenith. 

The number of the calendar month (universal time). 

The number of the day of the month including universal time as 
a decimal fraction. 

Local sidereal time. 

Weighting factor for the zenith angle Z to the meteor radiant. 

A constant factor in the exponent of the exponential weighting 
factor f 

d 

Weighting function for cumulative distribution of celestial 

latitude of radiant ~f f, f f ,. 

abed 

Data point for positive j3 . 

Data point for negative 

Weighting factor for obscured celestial latitude. 

Spatial weighting function ~ f^ f f^ f^ f^. 

Terrestrial weighting function ~ f f f , f„. 

& a c d f 

Unity, f , or f , depending on the analysis indicated, 
s t 

Average meteoroid material density in grams per cubic 
centimeter for photographic meteors and for smaller meteoroids 
with log m s -7.28. 

,X 5 Representation of log m , log v , /3 L e, and A, respectively. 

El I 6 1 

s. Sample weighted mean, deviation from the weighted mean, and 
1 weighted standard deviation for the variable where i — f , . . . ,5. 



DEFINITION OF SYMBOLS (Cont'd) 


Symbol 

r.. 

D 

R 


R.. 

D 


s 

e 


r 

1 • 2345 

Y 

ij . klm Partial correlation coefficient between X and X. with X^» X^.and 
X^ held fixed, where i, j, k, 1, m, = 1,! . . ,5. •* ' 

x Data point for log m not less than the weighted median log m. 

□ Data point for log m less than the weighted median log m. 

F> Mean number of sporadic and stream meteoroids per second 

per square meter of level surface with mass equal to or greater 
than m grams. 

■^M Mean number of sporadic and stream meteors per second per 

p square meter of level surface with absolute photographic magnitude 

equal to or less than M . 

P 

F Mean number of sporadic and stream meteoroids per second per 

square meter of level surface with momentum equal to or greater 
than mv gram kilometers per second. 

F 2 Mean number of sporadic and stream meteoroids per second per 
square meter of level surface with kinetic energy equal to or 
greater than mv 2 gram kilometers 2 second -2 . 


Definition 

Sample weighted correlation coefficient between X. and X. 
for i, j = 1 5. 1 J 

The determinant of the sample weighted correlation coefficients 

r.. for i, j = 1, . . . ,5. 
ij 

Cofactor of the element r.. in the determinant R. 

D 

Standard error of the estimate of X t from the equation of the 
regression plane. 

Multiple correlation coefficient for Xj in relation to X 2 ,. . . ,X 5 . 
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Symbol 

F mv ^ 2 


F, 

P 
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DEFINITION OF SYMBOLS ( Concluded) 

Definition 

Mean number of sporadic and stream meteoroids per second 
per square meter of revel surface with the geometric mean of 
momentum and kinetic energy equal to or greater than mv ^ 2 
gmkm ^ 2 sec-^ 2 - 

Factor by which the mean cumulative flux of sporadic and stream 
meteors exceeds the mean cumulative flux of sporadic mete ops. 

Negative of the radius exponent for the concentration of zodiacal 
light particles within a differential increment of radius . 

Designation of probable error. 
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OPIK'S LUMINOUS EFFICIENCY AND SUPPLEMENTED WHIPPLE WEIGHTING 

SUMMARY 


« • 

A formula for the mass of a meteoroid is derived from Opik's [ 1] physical 
theQry for dustball meteors, Whipple's [2] estimate of unit mass for a zero- 
absolute-visual-magnitude-thirty-kilometers-per-second meteor, and Hawkins 
and.Upton's [ 3 ] color index. After applying Dalton's [4] 0.20 order of magnitude 
decrement in the estimate of luminous efficiency due to higher pressures in the 
"U. S. Standard Atmosphere 1962" [5] and Opik's [6] formulation of known 
trail length, the masses and radiant celestial latitudes are calculated for Hawkins 
and Southworth's [7,8 ] random sample of 285 sporadic meteors. 

Weighting functions are derived for the statistical minimization of the 
effects of several sources of bias. Means, standard deviations, and ordinary, 
multiple, and partial correlations are shown for a 5-variable statistical analysis. 
Weighted cumulative distributions for parameters are plotted for two gradations 
with respect to mass. Data-point distributions are plotted which establish cumu- 
lative flux as functions of mass, momentum, kinetic energy, etc. , which show 
the effects of the earth's motion on the distributions and correlations of parameters. 

The results support the revisions in the author's [4] interpretation of 
Explorer XVI data and predicted puncture rates for the Pegasus meteoroid meas- 
urement satellites. The apparent discrepancy in the logarithmically linear 
relations between cumulative flux and mass, as reported by Whipple [ 2] for 
photographic meteor data and as reported by Beard [9, 10] for zodiacal light 
and solar F-corona data, is shown to be accountable through Opik's (see Whipple 
[ 11 ] ) probability that a particular meteoroid will encounter the earth in one 
revolution of the particle. 


SECTION I. INTRODUCTION 


A. JUSTIFICATION AND PURPOSE 

The technology of meteoroids and of their interaction with fields and 
with natural and artificial bodies in space is of considerable scientific and engi- 
neering interest. Each of the several sources of data continues to be of difficult 
interpretation due to indirection, extrapolation, and bias in random samples 



because the physical processes are selective. This report shows some new 
results for photographic meteors, from Opik's [ 1 ] physical theory for dustball 
meteors and from Hawkins and Southworth's [7,8] random sample of 285 sporadic 
meteors, which are in good agreement with Beard's [9,10] analysis of zodiacal 
light and solar F-corona data and with Dalton's [4] interpretation of Explorer XVI 
puncture data. 

B. METHOD AND NOTATION 

A multivariable statistical analysis is made with weighting functions 
of meteor height, velocity, celestial latitude, zenith angle, and earth-encounter 
probability. The sample is then equally divided with respect to an intermediate 
mass value, and weighted cumulative distributions for several parameters are 
plotted for the two gradations with respect to mass, jointly or separately. The 
analysis is repeated without the weighting with respect to the earth -encounter 
probability. All weighting factors are adjusted so that the sum for the sample is 
equal to the sample size. All flux values are cumulative with respect to the in- 
dicated parameter (e. g. , mass, momentum, etc.) and are in numbers per second 
per square meter of level surface. All logarithms are for base ten. Statistical 
notation is according to Hoel [ 12 ] . 


C. SCOPE 

The size of the data sample and the extent of analysis are not ex- 
haustive with respect to readily available data and techniques. Both could have 
been extended, but not conveniently, except with a larger computer than the GE 
225, which was used. Some known sources of bias remain uncorrected. No 
weighting is used with respect to seasonably and diumally varying sky brightness, 
with respect to magnitude above plate limit, nor with respect to either celestial 
longitude or longitude from the apex of the earth' s way , except insofar as it is 
related to the weighting with respect to zenith angle of radiant. 


SECTION II. PRELIMINARY ESTIMATE OF LUMINOUS EFFICIENCY 

FOR DUSTBALLS 


A. VISUAL AND PHOTOGRAPHIC MAGNITUDE 

Because of the nonlinear visual response to luminous intensity, sub- 
jective comparison of the luminous intensities of meteors and stars is facilitated 
by reckoning in magnitudes on a logarithmic scale such that an increment of five 
visual magnitudes decreases luminous intensity two orders of magnitude. The 


2 



relation between visual magnitude and photographic magnitude is complicated by 
the change in spectral sensitivity between cone vision and rod vision. Correspond- 
ing values given by Hawkins and Upton [ 3 ] , for the same ( Harvard) photographic 
meteor work from which the present data sample was taken, has been plotted in 
Figure i. When they are either very high or very low, the photographic and visual 
magnitudes can differ only by constants which depend on the kind of photographic 
emulsion which is used. Then, for the present data sample, 


M = M + 1.90, M ^ -2.25 

v P p 

= 0.725 M + 1.31, -2.25 < M < 1.52 
P P 


(1) 


= M + 0. 95, M s 1.52 
P P 

where is absolute visual magnitude and is absolute photographic magnitude. 


B. 


PHYSICAL THEORY OF METEORS 


as 


• • 

Opik 


1 1 ) gave a meteor reduction formula which can be expressed 


log m = k + log L - log 0 - 3 log v - ( 2/5) M , (2) 

V| V 

where k is a constant, L is the visible length of the trajectory, and 0 is the 
v l 

luminous efficiency. For dustballs, which he considered to constitute the major- 
ity of "ordinary" meteors, he considered the time of visibility t to be inversely 

proportional to the 0.7 power of the velocity v independently of meteoroid mass 
m ; i. e. , in equation ( 2) , 

log L = log v t v = k + 0.3 log v , (3) 

where k^ is a constant. Also in equation (2) , the meteor total luminous effi- 
ciency 0 for dustballs is the sum of the following two components: 

■ (4) 

where 0 q is the efficiency of meteor "initial" and "impace radiation," and 0^ is 

• • 

the efficiency of "temperature radiation. " Opik [ i ] shewed that for most visual 
meteors /? is the only important component. Therefore, log m in equation ( 2) 

can be approximated by 



where 


log m Q = k yi + l°g L “ log 0 Q - 3 log v - ( 2/5) . ( 6) 

• • 

Opik's [ 1 ] formula for the luminous efficiency main component /3 q is 

P Q = (Pi + Pl d p)/ (i + 9^) , (7) 

where / 3 t is the "initial fa' and "impact 0" efficiency for single meteor atoms 
decelerating in air without encounters with other meteor atoms, fa is the efficiency 

for secondary collisions between meteor atoms, and d. is the degree of dilution 

• • ' 

of the meteor coma. Opik [ 1] used one set of values of v in a table of values for 
9^m = 9 Q /ni and a different set of values of v in a tabulation for fa and fa. By 

combining his tables by linear interpolation with respect to log v, the values 
given in the upper half of Table I are found. 


• • 

Previously, Opik [6] showed that meteor luminous efficiency /3 is inversely 
proportional to the cube of the atmospheric pressure at the point of disappearance 
of the meteor. Opik's [ 1 ] values relating to luminous efficency, and therefore 
also all of the values in Table I, presuppose the Rocket Panel 1952 [ 13] model 
atmosphere. Instead of adjusting the values relating to luminous efficiency, the 
procedure which will be followed here is to increase the constant term in equation 
(2) by 0.20, which is three times the logarithm of the ratio of the means of the 
85-90 kilometer pressure from the U. S. Standard Atmosphere 1962 [5] and the 
ARDC Model Atmosphere 1959 [ 14 ] . The necessary adjustment due to the ARDC 
Model Atmosphere 1959 [14] is implicit in Whipple's [2] value for the meteoroid 
mass for the zero magnitude meteor which will be used to evaluate the constant 
term in equation ( 2) . 

The constant k in the equations ( 2) and ( 6) reduction formulas for visual 

meteors presupposes scotopic visual sensitivity and is therefore invariant only 
for limited values of visual magnitude (as Opik [ 1 ] cautioned) . For photographic 
meteors equations ( 2) , ( 3) , ( 5) , and ( 6) can be replaced by equations ( 8) and 
(9), 


logm =k -(2/5)M - 2.7 log v - log 0 (8) 

o Pi P ° 

log m = log m - log ( i + /3, //3 ) + k , (9) 

o t o p 2 
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TABLE I. 

VALUES OF PARAMETERS FOR CONSTRUCTION OF A METEOR 
REDUCTION FORMULA FOR DUSTBALLS 



log m 0 

= 3.81 log(Hj3 /0o) 




















































where k and k are constants which can be found from Whipple's [2 ] estimate 
Pi P 2 

that a i-gram meteoroid of velocity v = 30 corresponds to a meteor of zero 
absolute visual magnitude M^. With = -1.81 from equation ( 1) , with the 

values from the upper half of Table I at v = 30, and with equations ( 7) through ( 9) 
one finds k p by tentatively disregarding the distinction between and m; i.e. , 


k =0.106. (10) 

Pi 


The value of k is found by restoring the distinction between m and m; i. e. , 
P 2 0 


k = log ( 1 +P./P) @ (M , v)= (-1.81, 30). 

p 2 to p 


( 11 ) 


The numerical values in the upper part of Table I can be used to solve 
equation ( 8) for m Q at the indicated values of v and M . For this purpose , pa- 
rameters A^, Bg, and are introduced as follows: 


m 0 

B /3 " V m o 


( 12 ) 
(13) 


log = k - ( 2/5) M p - 2.7 log v - log ft . ( 14) 

Then by equations ( 7) and ( 12) through ( 14) , equation ( 8) can be solved explicitly 

for m as follows: 
o 

m o = {- (1 - B ^V + [ (1 - B /sV !+4A /3 C J‘ /2 } /2A g ' <15> 

The logarithms of the values of m^ from equation ( 15) are compiled in the 

lower part of Table I; and for each tabulated value, the logarithm of the corre- 
sponding luminous efficiency /3 , which by equations ( 7) , ( 8) , and ( 14) can be 
expressed by 


log j3 Q = - log + log m Q - log , ( 16) 

is plotted in Figure 2 to see how meteor luminous efficiency depends on velocity 

v. The smooth curve which seems most appropriate to represent the discrete 
points for -log /? is empirically chosen as 

G q (v) = 3.244 + 0.23 sin 5( log v - 1.55) . (17) 

The two points which fall somewhat below the given curve at the higher velocities 
are for the brightest meteors. 
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The discrepancies between the discrete values of log p in equation ( 16) 

and the corresponding functional representation in equation ( 17) can be incorpo- 
rated into the term which involves 6 in equation ( 9) as follows: 

G t ( M . v) = -log - log ( 1 + P^/PJ ~ G 0 (' v) ♦ ( 18) 

The discrete values for log ( 1 + PJj 3 ) which are shown in the lower part of 
Table I were computed from the other values in the same table by following 
Opik's [ 1 ] theory for the luminous efficiency of meteor "temperature radiation." 
The discrete values for the left side of equation ( 18) are shown in Figure 3 to be 
represented resonably well by the following linear least squares approximation: 

G ( M , v) = (1 - 0.14 M ) log v + 0.25 M - 1.534 , (19) 

t p ' P P 

• • 

within a maximum discrepancy of about 0.22 from Opik's [ 1] theory. 


By equation ( 11) the value of k from Table I would be 0.114, but with 

P2 

the further discrepancy of 0.020 from the least squares fit in Figure 3 at (M , 
= (-1.81 , 30) , one finds P 


v ) 
a' 


k = 0.134 . (20) 

P2 

Therefore, by equations ( 8) through ( 10) and ( 17) through ( 20) , one finds the 
reduction formula for photographic meteor data for ordinary (dustball) meteors 
based on the ARDC Model Atmosphere 1959 [ 14] ; 


log m = 1.95 - 1.7 log v + 0.23 sin 5( log v - 1.55) - ( 0.15 + 0.14 log v) M . ( 21) 

P 

C. ADJUSTMENT DUE TO U. S. STANDARD ATMOSPHERE 1962 

It was shown in Section II. B above that, by substituting the U. S. 
Standard Atmosphere 1962 [ 5] for the ARDC Model Atmosphere 1959 [ 14] , the 
estimated mass of a meteoroid which produces a given meteor must be increased 
by 0.20 order of magnitude. Therefore, equation (21) must be replaced by 

log m = 2.15 - 1.7 log v + 0.23 sin 5( log v - 1.55) - ( 0.15 + 0.14 log v) M . ( 22) 

P 
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D. VELOCITY AND LUMINOUS EFFICIENCY 

The construction in Figure 2 gives G q in equation (17) as an ap- 
proximation to -log /3 in equation (4) . Also, by equations (4) and ( 18) , the 
constructions in Figure 3 give G in equation ( 19) as an approximation to -log 

l 

/ 3 -G . Then, except for any constant component of log 0 which may have been 
o 

included in the constants k and k in equations ( 8) through ( 11) , and except 

Pi P2 

for the 0.20 constant decrement in log 0 due to a revision of the estimate of 
atmospheric pressure at typical meteor altitude ( see Section II. B) , log 0 is ap- 
proximated by - (G + G ) . Therefore, by equations ( 17) and ( 19) , the partial 
o t 

derivative of log 0 with respect to log v is 
n = 9 log 0/9 log v 

= - 1.15 cos 5( log v - 1.55) - (1.00 - 0.14 M ) . (23) 

• # 

Values from equation ( 23) are illustrated graphically in Figure 4. Opik 
[ 1 ] suggested further that the luminous efficiency of dustballs less than about 
one gram is approximately inversely proportional to velocity. This approximation 
would correspond to a unit negative value for n in equation ( 23) . The results in 
Figure 4 would seem to indicate that such a value for n would be rather typically 
a weighted average value for the present data sample. However, that approxi- 
mation is not used in the present analysis. 


SECTION III. METEOROID MASS WITH KNOWN TRAIL LENGTH 


The initial mass of a meteoroid, by equation (2) , is proportional to the 
detectable trail length L of the meteor which is produced in the atmosphere. The 
detectable trail length L, by equation ( 3) is equal to the product of the velocity 
v and the time of visibility t and is proportional to the 0.3 power of the velocity v. 

This result, which was used in the derivation of the equation (22) meteor reduction 
formula, was based on the following relation which Opik [ 1 ] gave for statistical 
use with dustball-type objects: 

t = 0.67 (21/v) °- 7 . (24) 
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But as a function of velocity v, zenith angle of meteor radiant Z, and the meteor 
heights at appearance and disappearance h^ , t is 

t =(h T1 -h )/vcos Z. (25) 

v B e ' ' 

• • 

Opik [ 6 ] previously stated that formulas for the masses of naked-eye meteors 
based on equation ( 25) are superior to those based on equation ( 24) except when 
the length of path or duration are unknown. By equations ( 2) and ( 3) it follows 
that equation (22) can be converted into the preferred meteor reduction formula 
for dustballs by subtracting the logarithm of equation ( 24) from the logarithm of 
equation ( 25) and adding the resulting null equation to equation ( 22) , giving 

logm = i.40 + log(h - h ) - log cos Z - 2 log v + 0.23 sin 5(log v - 1.55) 

-B 0 

-( 0.15 + 0.14 log v) M . (26) 

P 

Values for log m by equation ( 26) are tabulated in the second column of Appendix 
II. 


SECTION IV. ANALYSIS OF SPORADIC METEOR DATA 


A. DESCRIPTION OF THE DATA SAMPLE 

In 1958 Hawkins and Southworth [7 ] published data on a random 
sample of 360 of the meteors doubly photographed by the Baker Super-Schmidt 
cameras to a limiting photographic magnitude of +4 from stations at Donh Ana 
and Soledad, New Mexico, from February 1952 to July 1954. The sample was 
compiled according to the principle of the classical decimation process ( selecting 
every tenth detected meteor) . Of the 286 meteors which were judged to be 
sporadic, 285 were repeated by Hawkins and Southworth [ 8 ] in 1961 with a com- 
pilation of orbital elements and other data. These 285 meteors constitute the 
present data sample, with some values from each of references 7 and 8 , as 
tabulated in Appendix I . 

The two meteors of lowest velocity in Appendix I (lines 144 and 152) 
did not have data in reference 8 for elongation A and Whipple's cosmic weight 
f g . But for the other four meteors in Appendix I with velocity v s 11.9 (lines 

5, 61, 87, and 127) the average elongation and cosmic weight are 97.8 and 0.87 
respectively. These average values are assumed for lines 144 and 152. There 
are seventeen meteors in Appendix I (lines 12,14,20,21,159,178,185,192,195, 
197,205,206,209, 258, 280, 282, and 284) which did not have data in reference 8 for 
Whipple's cosmic weight f ; therefore, a unit value has been adopted in each case. 
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B. WEIGHTING FUNCTIONS 

Appendix I is thought to constitute a random sample of sporadic 
meteor data, but is also thought to have several components of bias due to physical 
selection effects , and therefore is not considered to represent a random sample 
of meteoroids in space until special weighting is considered. 

For data similar to the present sample, Whipple [ 11 ] used weighting 
factors ("cosmic weight") f proportional to P -1 v -2 for velocity v > 19 or pro- 
portional to 3 _1 P -1 v -2 for v < 19, where P is Opik's [ 15 ] probability that a 
meteoroid is a given orbit will encounter the earth in one revolution of the particle. 
McCrosky and Posen [ 16 ] and Jacchia and Whipple [ 17 ) used Whipple's [ 11 ] 
cosmic weight ~P -1 v -2 for all meteors; and McCrosky and Posen [ 16 ] said : 

"The uncertainty in the velocity mass law and in the number- luminosity law are 
such that attempted corrections for these effects will probably be in error by at 
least 1 in the velocity exponent. " 

It seems appropriate to consider that, for the population of meteoroids 
which encounter the earth, mass (m) may be statistically dependent on the 
velocity (v) , on the eccentricity (e) of the meteoroid heliocentric orbit, on the 
elongation (A.) of the true radiant from the apex of the earth's way, and on the 
celestial latitude (/3 ) of the corrected radiant. Also, the data is thought to be 
biased with respect to all of those parameters. The weighting factors which 
are developed in Sections IV. C through G below are thought to constitute a mini- 
mum necessary attention toward the reduction of bias for a significant statistical 
analysis. 

C. WEIGHTING FACTOR FOR HEIGHT AND VELOCITY 

Whipple's [ 11 ] cosmic weight has a factor v -2 which contains the 
factor v -1 / 2 , because the height was assumed to vary as v 1 / 4 and the effective 
collecting area was assumed to vary as h 2 . Since both velocity v and height h 
data are available, it is desirable to represent this weighting factor by 

f ~ h“ 2 v~ 3 / 2 . ( 27) 

3 . 

Values for this weighting factor, with sum equal to the sample size, are tabulated 
in the fourth column of Appendix II. 
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D. WEIGHTING FACTOR FOR EARTH-ENCOUNTER PROBABILITY 

Whipple's [ 11 ] cosmic weight f , which was taken from reference 

6 

8 and tabulated in Appendix I, is proportional to v -2 P -1 . Then a weighting factor 
f^ inversely proportional to Opik's [ 15] earth-encounter probability P is repre*- 

sented by 

f, ~ v 2 f ( 28) 

be 

Values for f^, with sum equal to the sample size, are tabulated in the fifth column 
of Appendix II. 

E. WEIGHTING FACTOR FOR CELESTIAL LATITUDE OF RADIANT 

Celestial latitude /3 is determined from the right ascension a 
e 

and the declination 6 (both tabulated in Appendix I) by the following formula 
from reference 18: 

/ 3 = sin -1 ( 0.91741 sin 6 - 0.39795 cos 6 sin a) . (29) 

6 

Values for /? by equation ( 29) are tabulated in the third column of Appendix II. 

It seems appropriate that meteors which have radiants on circles of 

celestial latitude with not more than 95 percent occlusion (by the horizon) should 

be weighted by a factor f in near inverse proportion to the apparent ( above the 

c 

horizon) fraction. The weighting factor for reducing bias with respect to the 

celestial latitude of the radiant is therefore related to the difference 2 A between 

e 

the celestial logitudes of the horizon points on the circle of celestial latitude fi 

6 

through the meteor radiant. 

The pole of the ecliptic, the zenith, and either one of the horizon points 

on the circle of celestial latitude /3 forms a spherical triangle such that the angle 

© 

at the pole of the ecliptic is tt - X with opposite side ?r/2’and with adiacent sides 

e 

7r/2 - /3 and 7 r/2 - (3 where j3 is the celestial latitude of the zenith. The solution 
e z z 

to the spnencai ixinu^ic v cs 

X =0 when 

e 

cos X = tan j8 tan 0 when 

e e z 


P e 2 n/2 - 

< »/ 2 - e z 


<30) 
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At the Dona Ana station tan (3 is always positive, 0 ^ A — tt/2 when (3 >0, 

Z 6 G 

7r/2 £ A ^ 7 T when /3 < 0, and the circles of celestial latitude (3 > 7r/2 - (3 are 
e e e z 

entirely above the horizon. Therefore, 

f = tt/( 7r - X ) when /3 < ir/2 - /3 and n/ (it - X ) ^20 

c e' e z e' 

= 20 when /3 < 7r/2 — /3 and tt/ (it - X )> 20 

G Z G 

Kells, Kern, and Bland [19] state that the right ascension a of the 

z 

zenith of a place is the sidereal time at that place; and Russell, Dugan and 

Stewart [ 20] state that sidereal time is usually correct within four or five minutes 

when approximated from civil time by assuming that on September 22 the two 

times agree and that sidereal time gains two hours each month and proportionally 

each day and hour. Also, Hawkins and Southworth [ 7] state that to obtain the 

time of appearance in local mean time at the Dona Ana station, 0.29667 must be 

subtracted from the date given ( in Appendix I ) ; and Hawkins [ 21 ] gave the . 

declination 6 of the zenith of that station as 
z 

6 = sin -1 0.8433343 = cos -1 0.5373893 . (32) 

z 

With this information, together with equation ( 29), the celestial latitude (3 of 

Z 

the zenith can be formulated. Local sidereal day t , at the Dona Ana station is 

ds 

approximated by 

t =t J + (t - 9) ( 2/24) + ( t . - 22)/ 365.26 - 0.29667 (366.26/365.26) , 
ds du mu du. (33) 

where t and t , are the numbers of the month and day universal time, respec- 
mu du 

tively ( Appendix I) . The right ascension a of the zenith is given in radians by 

z 

ct = 2tt t , ( 34) 

z s 

where t is the local sidereal time expressed as a decimal fraction of a day; i. e. , 

t is the decimal residue which is left after subtracting the next lower integral 
s 

value (positive, nil, or negative) from the local sidereal day t^ , equation (33). 

By equations (29) , (32) , and ( 34) , the celestial latitude (3 of the zenith is re- 

z 

lated to the local sidereal time t by 

s 
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sin /3 = 0.77368 - 0.21385 sin 2 tt t 
z ! 

34.0 ^80.9° 

z 

Equations (33) through (35) support the solution of equations (30) and 

( 31) for the weighting factor f with respect to meteor radiant celestial latitude. 

c 

The results, with values normalized so that their sum is equal to the sample 
weight, are tabulated in the sixth column of Appendix II. 

F. WEIGHTING FACTOR FOR ZENITH ANGLE OF RADIANT 

It has not been established whether or not there may be some bias 
directly for the zenith angle to the meteor radiant. Elford, Hawkins, and 
Southworth's [22] analysis of radio meteors indicates that, after corrections 
for antenna selectivity and velocity selection effects, radiants are markedly non- 
uniformly distributed in elongation A from the apex of the earth’s way. Indirectly, 
this effect will bias the distribution of zenith angles Z. 

Part of the bias in elongation A can be reduced by weighting with respect 
to zenith angle Z, as follows: 
k Z 

f d = e Z , (36) 

where k is a constant to be determined from the distribution of meteor celestial 
z 

latitude /3 . 

e 

A first approximation to k in equation ( 36) is made by noting that the 

Z 

uniformly weighted sample mean zenith angle Z is 0.64 instead of the value 7r/4 
which would be expected from an isotropic distribution. Then, 

-k Z -k Z 

0.64= f n/2 e Z Z sin 2ZdZ/ f n/2 e Z sin 2ZdZ 
J o *o 

-k tt/2 f -k tt/2\ 

= (tt/2) e Z Ai+e Z )+ 2k/(k 2 + 4) (37) 

z z 

= 2.972 - il. 707 k + 6.053 k 2 + . . . 

z z 

The solutions to the linear and quadratic truncations of the Maclaurin series in 
equation (37) are 0.20 and 0.225, respectively. But when 0.225 is substituted 
for k^ in equation ( 36) , and the weighting function f^ which is used in finding the 

sample cumulative distribution of /3 is 
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( 38 ) 


f ~ f f f f . , 

o abed 

the results illustrated in Figure 5 are found. Figure 6 shows that a more nearly 
symmetric distribution of /3 results for vanishing k . Also Figure 7 shows the 

6 Z 

results for k = -0.18, which is the value which minimizes the sum of the squares 
z 

of the differences between the cumulative distributions for positive $ and negative 

G 

B for -42° ^ /3 £ 42°. The latter results seem more reasonable; and the value 
e e 

k = -0.18 is adopted for the rest of the analysis. The corresponding values of 
z 

f and f in equations ( 36) and ( 38) , normalized so that their sum is equal to the 
d o 

sample size, are tabulated in the seventh and eighth columns of Appendix II, 
respectively. 


G. 


WEIGHTING FACTOR FOR OBSCURED CELESTIAL LATITUDE 


In the rest of this analysis, the celestial latitude of the meteor 
radiant will be treated only arithmetically. Because the observations were made 
in New Mexico, radiants with celestial latitudes between -42 and -90 degrees 
were not subject to observation. By inspection of Figure 7, it seems reasonable 
to minimize the bias in the cumulative distribution of |/3 fi j by doubling the weight 
for the meteors with celestial latitudes between 42 and 90 degrees. The purpose 
is served with the following weighting factor: 


f f = 1 for |/3 e | ^ 42 


= 2 for p >42 
e 


H. WEIGHTED STATISTICAL ANALYSIS 


(39) 


The results of a multivariable statistical analysis, performed 

three times with different weighting functions, are tabulated in Appendix HI. 

Each line in the tabulation is for a different statistical parameter, identified in 

the first column. The results for uniform weighting are given in the second 

column. The values obtained with weighting functions f and f. are given in the 

s t 

third and fourth columns, respectively, where 
f ~f f, > 

S ° f l (40 

~f £f f.f. ( 

a b c d f J 
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and 


f /f. 
s b 

'f f f. 
a c d f 


( 41 ) 


The function f is for an observer with heliocentric orbital motion similar to 

V 

earth's, and therefore does not involve the meteoroid-earth encounter probability. 

The function f is for considerations of bodies in space , whether or not they may 
s 

encounter the moving earth. The values of f and f in equations ( 40) and ( 41) , 

s t 

no localized so that their sum is equal to the sample size, are tabulated in the 
ninth and tenth columns of Appendix II, respectively. 


The indicated mass m for some of the meteors in Appendix n may be 
somewhat in error because of the assumption that all of the meteors are due to 
dustballs. Before proceeding with the statistical analysis, it is appropriate to 
estimate the relative number of meteors in the sample which may not be due to 
dustballs. 


Opik [ 1 ] gives the material density of meteoritic stone as 3.4 grams per 
cubic centimeter. Whipple and Hawkins [ 23 ] found that photographic meteors 
are produced by bodies with material density approximately 0.1 with a maximum 
value of 0.3; and they concluded ( 1) that radio and visual meteors are almost 
entirely of cometary origin; ( 2) that detonating bolides and meteorites are per- 
haps entirely of asteroidal origin; ( 3) that "... if asteroidal meteors are in fact 
present among the sample of photographic meteors, they number at most a very 
few percent, probably fewer than 3 percent;" and (4) that the low value of mate- 
rial density "... is not surprising in view of the fragile nature of meteoroids as 
exhibited by the phenomenon of fragmentation. " Later, in their analysis of pre- 
cision orbits of photographic meteors, Jacchia and Whipple [ 17] stated: "The 
writers are of the opinion that the asteroidal contribution to the photographic 
meteors probably does not exceed 1 percent of the total and may well be less. " 

Later yet, Briggs [ 24] found that an average material density of p of 

P 

o p = 0.1 , (42) 

(for log m ^ -7.28 i. e, for particles with radii not less than 50 microns) , is con- 
sistent with the interpretation of zodiacal-light observations and the theoretical 
apparent brightness due to scattered sunlight from the steady-stale system of 
particles under Poynting- Robertson action. Therefore, with such a low average 
density, it does not seem likely that a random sample of 285 sporadic meteors 
would include many stones or irons, if any. 
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Except for the use of weighting factors, the mathematical treatment in 
this section, including terminology and symbols, is according to Hoel [ 12 ] . The 
weighting function, unity, f g , or f , as the c.ase may be, will be designated by 
f in the formulation below. 

Let the 5 capital letters, Xj,. . . ,X 5 designate the variables for which a 
statistical analysis is desired, as follows: 

Xj, . . . ,X 5 = log m , logv a ,|^ e j , e, and X, respectively. (43) 

The sample weighted means X., deviation x. from the sample weighted means, 
and the sample weighted standard deviations s. for the variables X^ are given by 


= (1/285)2 fX. , 

(44) 

= X. -X. , 

l l 

(45) 

= [( 1/285) 2 fx\ ] ^ , 

(46) 


where i = 1 , . . . , 5. 


The sample weighted correlation coefficients r.. between the variables X. and 
X^ are given by ^ 

r =r = ( 1/285 s.s.) 2 fx.x. where i, i = 1,. . . ,5. (47) 

i] Ji i] i 3 


The determinant R of the sample weighted correlation coefficients r.. is 
I r ll • • • r 15 ' 


R = 


r 5 i 


r 55 


(48) 


The equation of the least squares regression plane for the purpose of estimating 
X t is 


Rn (xj/si) + ... + R 15 (x 5 /s 5 )= 0 , 


(49) 
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where R.. is the cofactor of the element r.. in R; i.e. , R.. is (-1) 1+ J times the 
ij ij ij 

minor of r.. inR. The standard error of estimate s , where X, is estimated 
ij e 

from the equation of the regression plane, and the multiple correlation coefficient 

2345 are 

s = Sl (R/Rn) (50) 

V 


and 


r l* 2345 - C 1 “ ( R/Hll)] 


V, 


(51) 


The partial correlation coefficient r.. . . between X. and X. with X, , X, , and 

ij • klm i j k 1 

X held fixed is 
m 

r _ _ r / g 1/2 

ij • klm ij ' ii jj ; where i, j, k, 1, m= 1 ,. . . ,5. (52) 

The results of this 5-variable statistical analysis are shown in Appendix 
in, based on the we ighting factors from Appendix II. 


I. 


SAMPLE PARTITION FOR TWO REGIMES OF MASS 


The meteor for line No. 109 in Appendix II is found to have very 
nearly the weighted median value of log m. All of the sample meteoroids which 
do not have smaller mass are considered "large" and are given crossed data 
points in Figures 10 through 21. The other meteoroids are said to be "small" 
and are given square data points in Figures 10 through 21. In Figures 11 through 
15, the weighted sample cumulative distributions of parameters are shown for 
these two regimes of mass separately. 

J. WEIGHTED SAMPLE CUMULATIVE DISTRIBUTIONS OF 
PARAMETERS 


Scatter diagrams, such as Figures 8 and 9, help to show how a 
pair of parameters are related; e. g. , Figure 8 shows to what extent faster 
meteors tend to be brighter, and Figure 9 shows how the faster ones tend to have 
radiants closer to the celestial equator. But the relations between three param- 
eters, as in Figure 10, are more difficult to interpret visually, even with no 
consideration to nonuniform weighting. The weighted relation between pairs of 
parameters is shown more clearly, as in Figures 11 through 15, by the weighted 
sample cumulative distributions for one parameter over separate regimes ( strata) 
of the other parameter. 
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Figure ii shows the "high mass" and "low mass" separate cumulative 

distributions of orbital eccentricity with "spatial" weighting f . The results apply 

s 

to meteoroids in interplanetary space, regardless of whether or not they may 

encounter the earth. Figure 12 is similar to Figure 11 except that f is replaced 

s 

by "terrestrial" weighting f ; but the results apply to those meteoroids which 

l 

actually encounter the earth, or encounter an observer in an orbit similar to the 
earth's. Figures 13 and 14 show for the celestial latitude of the meteor radiant 
what Figures 11 and 12, respectively, show for eccentricity. Figure 15 shows 
for the air-entry velocity what Figures 12 and 14 showed for the eccentricity and 
the celestial latitude of the radiant, respectively. 

Weighted sample cumulative distributions for several parameters, indi- 
vidually, are of considerable interest because, when the surveillance area and 
duration are known for the data sample, they permit the determination of weighted 
mean cumulative flux with respect to those parameters. When the weighted sample 
cumulative distribution is plotted as a function of the parameter of interest, the 
slope is the same as the slope of the weighted mean flux; and the slope of the 
logarithm of the same cumulative distribution is the same as the slope of the 
logarithm of the flux. Figures 16 through 21 show automatic plots of the logarithm 
of the sample cumulative distribution of the two strata of the sample collectively, 
with different weighting and for different observational or dynamic parameters, 
or functions of parameters possibly of direct interest for models of particle con- 
centration or of impact and puncture fluxes. 

Because the sample size is 285 and is equal to the sample sum of the 
values of the weighting function, the logarithm of the weighted cumulative distri- 
bution is log (285/2) = 2.15 at the second quartile. For each case in Figures 16 
through 21 the sample weighted cumulative distribution is seen to be approximately 
logarithmically linear up to or somewhat beyond the second quartile. In choosing 
a line of best fit, the upper half of the sample is ignored because the weighting 
functions do not include any factor with respect to magnitude-above-plate limit. 

K. PRELIMINARY ESTIMATES FOR MEAN CUMULATIVE FLUX 
OF SPORADIC METEORS 

In 1957, with similar but somewhat smaller random sample of 
sporadic meteors ( 247 instead of 285) , Hawkins and Upton [ 3, 25] found that 
the logarithm of the weighted mean number of sporadic meteors per square kilo- 
meter per hour with absolute photographic magnitude equal to or less than 

is equal to 0,537 M - 4. 33. Figure 18 shows that a line with a 0.537 slope agrees 
P 
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well with the present results. Therefore, Hawkins and Upton's [3 ] results for 

flux versus M are used to infer from the results in Figure 18 that the area-time 
P 


product for the present data sample is 10 15 * 18 square meter seconds (i. e. , 5980 
square kilometers with plate-exposure intervals totaling 70.2 hours) . 


The logarithm of the sample cumulative distribution is approximated by 
the logarithm of the population mean cumulative flux plus the logarithm of the 
area-time exposure product. Therefore, the lines fitted in Figures 17 through 
21 imply the following equations (53) through (57) respectively, for the mean 
cumulative sporadic meteor flux per second onto one square meter of level sur- 
face at the top of the atmosphere: 


log F> - log F = -1.34 logm - 14.52 

(53) 

log F„ - log F = 0.537 M - 13.89 
M p 

P 

(54) 

log F - log F = -1.09 log (mv) - 12.86 

mv 

(55) 

log F mv 2 - log F = -0.92 log (mv 2 ) - 14.00 

(56) 

log F mv 3 / 2 " l0g F = " l0S (mvV2) " i3 * 92 ’ 

(57) 


where F is the factor by which the mean cumulative flux of sporadic and stream 
meteors exceeds the mean cumulative flux of sporadic meteors. 


SECTION V. MEAN FLUX COMPONENT FOR STREAM METEORS 


Hawkins and Southworth [ 7 ] published a random sample of 74 stream 
meteors along with their random sample of 286 sporadic meters. In referring 
to this work of Hawkins and Southworth [ 7 ] , Whipple and Hawkins [ 23 ] indicated 
that 83 percent of such meteors do not belong to major streams. Therefore, a 
model for the mean total flux would seem to be obtained by increasing the mean 
flux for sporadic meteors by 0.08 order of magnitude; i. e. , in equations (53) 
through ( 57) , 

log F = -log 0.83 = 0.08 . ( 58) 
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SECTION VI. REVISED ESTIMATE OF FLUX FROM EXPLORER XVI 

PUNCTURE DATA 


Dalton's [4] analysis of Explorer XVI puncture data indicated that the 
flux of meteoroids with mass equal to or greater than m grams onto a randomly 
oriented surface just above the atmosphere would be exp 10 (-i.34 logm - 14.92) . 
T his result involved a puncture flux enhancement factor, 4 = iO 0,60 . But in a 
more rigorous later derivation, Dalton [ 26] showed that the puncture flux en- 
hancement factor is 5 = i O 0 * 70 ; i.e. , the incident flux should have been estimated 
0. 10 order of magnitude higher. Also, the mean value of log v with spatial 

weighting f which was available for those efforts is higher than the correspond- 
s 

ing value with terrestrial weighting f ; i.e. , from Appendix HI the two weighted 

l 

means of X 2 = log v arelc 2 = log 26.7 and log 19.4, respectively. And, according 
to the puncture flux model which was used for those efforts, mass sufficient to 
puncture a given structure varies inversely with the 3/2 power of the closing 
velocity, and incident flux varies inversely with the 1.34 power of the nominally' 
sufficient mass. Therefore, the estimated nominally sufficient mass for punc- 
turing the Explorer XVI transducers should be increased by 0.21 order of mag- 
nitude, and the estimated incident flux should be increased by 0.28 order of mag- 
nitudes due to the lower estimated velocities. Then, when the 0.10 and 0.28- 
order-of-magnitudes adjustments are combined with the 0. 30-order-of-magnitude 
increase in flux due to removing the earth shielding factor, the estimate for log 
F> which is supported by the Explorer XVI results becomes exp 10 (-1. 34 logm -14. 24) . 


SECTION VII. REVISED ESTIMATES OF MASS AND FLUX FROM 
PHOTOGRAPHIC METEOR DATA 


The revised estimate (in Section VI above) for log F > based on Explorer 

XVI puncture data is 0.20 order of magnitude higher than the result by adding 
equations ( 53) and ( 58) . These results can be brought into agreement by in- 
creasing the estimated mass of the photographic meteors by 0.15 order of mag- 
nitude. The mass of the meteoroid which produces a zero absolute visual mag- 
nitude meteor by entering the air at 30 kilometers per second would therefore 
be 2.24 = lo 0, 35 instead of 10°* 20 grams. Although this adjustment is well within 
the uncertainty of the results, it is supported by the value 4.4 grams which 
Hawkins [ 27] reported in 1964, By equations ( 53) and ( 55) through ( 57) , the 

order-of-magnitude increases in F , F , F 2 , and F o/ 2 , due to this 0.15- 
6 > mv mv^ mv J/ 

order-of-magnitude increase in the estimate of meteoroid mass, are 0.15 times 

1,34, 1.09, 0.92, and 1.00, respectively ; i.e. , with equations (58) , 
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(59) 


log F = - 1*34 log m - 14.24 ± 0.60 

log F,. = 0.537 M - 13.81 ± 0.15 (60) 

& M p 

P 

log F^^ = - 1,09 log (mv) - 12.62 ±0,60 (61) 

log = - 0.92 log (mv 2 ) - 13.78±0.60 (62) 

!og F 3 / 2 = - log (mv 3 / 2 ) - 13.69 ±0.60 (63) 


SECTION VIII. STATISTICAL SIGNIFICANCE 
AND PHYSICAL IMPLICATIONS 


The direct correlation between log m and log v, with terrestrial weighting 
f^, is only 0.01, as tabulated in Appendix III. But the corresponding partial 
correlation is 0.10. Although numerically small, this partial correlation is 
quite significant in changing the slopes of the log cumulative distributions with 
respect to mass, momentum, kinetic energy, or other functions of mass and 
velocity, as shown in Figures 17, and 19 through 21. 

The slopes of the log cumulative distributions with respect to momentum, 
kinetic energy, etc. , are established just as accurately in Figures 19 through 21 
as is established in Figure 17 with respect to mass. Therefore, equations (59) 
through (63) , which are based on these results, establish more substantial models 
for extrapolation than one would have by assuming that mass and velocity were 
statistically independent. 

Figure 17 establishes through equation (59) that the flux of meteoroids, 
with mass not less thanm, which encounter the earth, varies inversely with 
the 1.34 power of m. But, by Figure 16, the corresponding concentration of 
meteoroids in space, which may or may not encounter the earth, varies inversely 
with the 1.00 power of m. This change in the mass exponent results by dividing 
the terrestrial weighting function f^ by Opik's [ 15] earth-encounter probability P. 

This result is significant in the physical interpretation of zodiacal- light data. 

The mass exponent for the concentration of particles with mass equal to 
or greater than m , times the radius exponent in the expression for mass, is 
equal to -p + 1, where -p is the radius exponent for the concentration of particles 
within a differential increment of radius. Briggs [24] indicates that particles in 


21 


the 1 to 50 micron range play the dominant role in the zodiacal light; and Ehricke 
[ 18] says that, according to Takakubo, the radius of the particles most effective 
is about 20 microns. Also Briggs [24] finds that the material density of particles 
with radii between 0.6 and 50 microns should be considered inversely proportional 
to the 0.678 power of the radius; i. e. , mass varies with the 3 -0.678 = 2. 322 
power of the radius. Therefore, the unit negative slope in Figure 16 corresponds 
to 


p = 3.32 . (64) 

In 1958 Beard [ 9 ] reported: "The analysis of zodiacal light . . . results in p > 3 
. . . possibly 3.5 ... "; and in 1963 Beard [ 10] indicated a preference for the 
value p = 3.5. Actually, there is an uncertainty in the line of best fit in Figure 16 
such that a line with -1.08 slope would seem to fit as well as the one which is 
illustrated with slope -1.00. This would correspond to p = 3.5 if material density 
varies, according to Briggs [24] , inversely with the 0.678 power of the radius. 
Also, the unit negative slope in Figure 16 and the value p = 3.5 could both be 
correct if material density varies inversely with the 0.500 power of the radius. 
Anyway, it is not necessary to assume that the mass dependence of particle con- 
centration in interplanetary space is different for zodiacal-light particles and 
photographic meteors with seven or eight orders of magnitude larger mass. 


SECTION IX. REVISED PREDICTION OF PUNCTURE RATES 
IN PEGASUS EXPERIMENTS 


In January 1965 Dalton [4] published predicted puncture rates for the 
1.5-mil aluminum 1100-H14 and the 8 and 16-mil aluminum 2024-T3 puncture 
transducers on the Pegasus meteoroid measurement satellites. The effect of 
the copper foil plate of the capacitor, thermal control coating, trilaminate of 
mylar, and layers of adhesive were considered to increase the puncture resistance 
equivalent to a 0.4-mil increase in the thickness of the aluminum foil or sheet. 
Puncture fluxes for Pegasus transducers were scaled from the 44 punctures in 
10 7,06 ° square meter seconds ( 1431 square foot days) , and 11 punctures in 10 6,738 
square meter seconds (682 square foot days) , for the Explorer XVI puncture 
transducers of 1.15 and 2.15-mil beryllium copper, respectively, thought to be 
equivalent to twice those thicknesses of aluminum 2024-T3. In that extrapolation, 
it was erroneously anticipated that the slope illustrated in Figure 21 would be the 
same as in Figure 17. Also, no consideration was given to the fact that the 
Explorer XVI puncture transducers are thin targets, unlike the Pegasus trans- 
ducers which are backed with a thick foam material and therefore may respond 
more nearly like semi-infinite targets. 
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According to the hypervelocity impact model which the author developed 
in reference 28, a back-supported sheet should have the same puncture resistance 
as a sheet 59 percent thicker and of the same material without a back support; 
and the cube of puncture depth should be proportional to mv 3 / 2 . Then, by equation 
(63) , the ratio of puncture fluxes through sheets of different thicknesses of the 
same material should be the cube of the reciprocal of the thickness ratio. Also, 
by Reismann, Donahue, and Bu^kitt's [29] multivariable analysis of hypervelocity 
impact data, a given sheet of aluminum 2024- T3 should have the same puncture 
resistance as a 4.3 percent thicker sheet of aluminum 1100-H14. Therefore, the 
puncture transducers on the Pegasus satellites are thought to have the same 
puncture resistance as the following thicknesses of aluminum 2024-T3 sheets 
without back- supporting: 2.9, 13.4 and 26.1 m-ils. 

If one assumes that the earth-shielding factor is approximately the same 
for Explorer XVI and Pegasus orbits, then the puncture fluxes for the three 
Pegasus sensors should be 10 -5,719 , IO -7713 , and io -8 * 582 per square meter per 
second, based only on the 44 punctures in the thinnest material on Explorer XVI. 
These fluxes should be io -5,184 , 1CT 7 * 178 , and io -8,046 per square meter per 
second, based only on the 11 punctures in the thicker material on Explorer XVI. 
Although 44 punctures are somewhat more statistically significant than 11, the 
predictions based on the 11 punctures involve less extrapolation for the two thicker 
transducers on Pegasus. Therefore, the geometric means will be assumed as 
IO -5,451 , io -7 * 445 , and io -8,314 per square meter per second. Because each panel 
on Pegasus has an area of 0.516 square meters, the predicted number of punctures 
per panel per day are 0.158, 0.00160, and 0.000216. Then, if all of the panels 
would remain responsive f&r one year, the punctures should be 920, 20, and 29 
for the 1.5, 8, and 16-mil panels, respectively. Sufficient flight data have not 
yet been evaluated (as of June 8, 1965) for checking the accuracy of this predic- 
tion. 


SECTION X. CONCLUSIONS 


The most pertinent results are summarized in Figure 22. More detailed 
results of the 5-variate statistical analysis, with the input data from Appendix I 
and the weighting factors from Appendix II. are tabulated in Appendix EH. In 
addition to the cumulative distributions which are plotted in Figures 16 through 21, 
supporting the flux formulas which are summarized in Figure 22, Figures 11 
through 15 show how the cumulative distributions of parameters depend on the 
stratification with respect to mass and on the inverse weighting with respect to 
the earth-encounter probability. 
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By discarding the earth-encounter probability ( substituting the terrestrial 

weighting factor f for the spatial weighting factor f in Appendix II) , mean log v 
t s 

is decreased from log 26.7 to log 19.4, mean |'/3 e | is decreased from 28.9 to 25.9 

degrees, mean eccentricity is decreased from 0.70 to 0.55, and the mass exponent 
for the concentration or flux of particles of equal or greater mass is changed 
from -1.00 to -1.34. 
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Mp, Meteor Absolute Photographic Magnitude 
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FIGURE i. PHOTOGRAPHIC SENSITIVITY AND THE PHOTOPIC, PURKINJE, AND 
SCOTOPIC VISUAL SENSITIVITIES FOR METEORS 2 
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FIGURE 5. CUMULATIVE DISTRIBUTION OF CELESTIAL LATITUDE OF 
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Weighted Sample Cumulative Distribution for 
either Positive or Negative /J e Separately 
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Weighted Sample Cumulative Distribution for 
either Positive or Negative /3 e Separately 



FIGURE 7. CUMULATIVE DISTRIBUTION OF CELESTIAL LATITUDE OF 

RADIANT WEIGHTED exp (-0. 18Z) FOR ZENITH ANGLE Z 

6 


31 



NITUDE 





h. Height at Maximum Brilliance 



FIGURE 10. METEOR HEIGHT VERSUS ZENITH ANGLE FOR TWO REGIMES 
OF MASS 
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FIGURE li. SPATIALLY WEIGHTED DISTRIBUTIONS OF ECCENTRICITY 
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Sample Cumulative Distributions Separately for High and Low 
Mass Regimes with Terrestrial Weighted f* 
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Sample Cumulative Distribution Separately for 'High' and 'Low' 
Mass Regimes with Spatial Weighting f s 



FIGURE 13. SPATIALLY WEIGHTED DISTRIBUTIONS OF CELESTIAL LATITUDE 
OF RADIANT 
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Sample Cumulative Distributions Separately for 'High' and 'Low' 
Mass Regimes with Terrestial Weighting ft 
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FIGURE 14. TERRESTRIALLY WEIGHTED DISTRIBUTIONS OF CELESTIAL 
LATITUDE OF RADIANT 


38 





Sample Cumulative Distributions Separately for 'High' and 'Low' 
Mass Regimes with Terrestrial Weighing ff 



FIGURE lt>. TERRESTRIALLY WEIGHTED DISTRIBUTIONS OF AIR-ENTRY 
VELOCITY 
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Logarithm of Sample Cumulative Distribution 



FIGURE 16. SPATIALLY WEIGHTED DISTRIBUTION OF METEOROID MASS 
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Logarithm of Sample Cumulative Distribution 
with Terrestrial Weighting ff 



FIGURE 17. TERRESTRIALLY WEIGHTED DISTRIBUTION OF METEOROID 
MASS 
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Logarithm of Sample Cumulative Distribution 
with Terrestrial Weighting tf 



FIGURE 18. TERRESTRIALLY WEIGHTED DISTRIBUTION OF ABSOLUTE 
PHOTOGRAPHIC MAGNITUDE 
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Logarithm of Sample Cumulative Distribution 
with Terrestrial Weighting ft 



FIGURE 19. TERRESTRIALLY WEIGHTED DISTRIBUTION OF METEOROID 
AIR- ENTRY MOMENTUM 
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Logarithm of Sample Cumulative Distribution 
with Terrestrial Weighting 



FIGURE 20. TERRESTRIALLY WEIGHTED METEOROID AIR-ENTRY KINETIC 
ENERGY 
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Logarithm of Sample Cumulative Distribution 
with Terrestrial Weighting f| 
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TERRESTB.IAT.TiY WEIGHTED DISTRIBUTION OF THE GEOMETRIC 
MEAN OF METEOROID AIR-ENTRY MOMENTUM AND KINETIC 
ENERGY 
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2.24 grams mass for zero absolute visual magnitude meteor 
at 30 kilometers per second 


log 

F Mp = 

0.537 

M p -13.8 

1 ± 0.15 

log 

F > = - 

1 .34 

log m - 14. 

24 ± 0.60 

log 

F mv = 

-1.09 

log (mv) - 1 

2.62 ± 0.60 

log 

F mv* = 

-0.92 log (mv*) 

-13.78 ± 0.60 

log 

Fmv 9 '* 

= - log (mv s/t ) - 

13.69 ± 0.60 


PARTIAL CORRELATIONS 

Terrestrial Weighting f^ Statial Weighting f s 

0.097 for log m vs. log v 0.107 for log m vs. |/9 e | 

0. 103 for log m vs. |£ e | -0.200 for log m vs. e 

-0.403 for log v vs. |£ e | 0 516 for l#e| vs - e 

0.882 for log v vs. X -0.799 for X vs. e 
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FIGURE 22. SUMMARY RESULTS 
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APPENDIX H. COMPUTED VALUES IN THE DATA SAMPLE 
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APPENDIX H. (Concluded) 






APPENDIX III. NUMERICAL RESULTS FROM THE MULTIVARIABLE 

STATISTICAL ANALYSIS 


X lf . 

. . ,X 5 = log m, log v . 

j /3 e j , e, and A, 

Respectively 

Statistical 

Parameter 

With Uniform 
Weighting 

With Spatial 
Weighting f 

s 

With Terrestrial 
Weighting f 

u 

x t 

-1.0351 

-1.0503 

-1.0649 

X 2 

1,4666 

1.4240 

1.2870 

Xg 

26,0074 

28.9383 

25.8878 


0,7242 

0.7040 

0.5505 

X 5 

78.1870 

83.1402 

97.1172 

S 1 

0.4748 

0.4668 

0.4170 

s 2 

0.2287 

0.1793 

0.1684 

s 3 

22.3558 

18.9195 

22,2727 

s 4 

0.2200 

0.2214 

0.2418 

s 5 

31.6016 

21.2353 

24,4780 

r ll 

1.0000 

1.0000 

1.0000 

r 12 

0.0241 

0.0551 

0.0102 

r 13 

0.0028 

0.0040 

0.0424 

r 14 

0.1122 

0.1972 

0.1148 

r 15 

0.0335 

0.0770 

0.1082 
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APPENDIX III, 


( Cont'd) 


Statistical 

Parameter 

With Uniform 
Weighting 

With Spatial 
Weighting f 

s 

With Terrestrial 
Weighting f 

l 

r 22 

1.0000 

1.0000 

1,0000 

*23 

-0.0485 

-0. 1242 

0.0465 

*24 

0.7917 

0.8324 

0.8108 

*25 

-0. 8963 

-0.8258 

-0.6369 

*33 

1.0000 

1.0000 

1.0000 

*34 

-0.1969 

-0.2718 

-0.1343 

*35 

0.0829 

0.1196 

-0.0606 

*44 

1.0000 

1.0000 

1.0000 

*45 

-0.5186 

-0.4496 

-0.1649 

*55 

1.0000 

1.0000 

1.0000 
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-0.0046 
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0.0075 

0.0189 

0.0273 
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APPENDIX in. (Cont’d) 


Statistical 

With Uniform 

With Spatial 

With Terrestrial 

Parameter 
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s 
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»ll 

-0.0006 
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R 22 
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APPENDIX III. (Concluded) 


Statistical 

Parameter 

With Uniform 
Weighting 

With Spatial 
Weighting f 
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With Terrestrial 
Weighting f 
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r 14* 235 

-0.1002 
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-0.1401 

r 15- 234 

0.0060 

-0.0321 
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r 23- 145 

-0.4300 
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0.4576 
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r35* 124 

0.3965 
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r 45> 123 

-0.7525 

-0.7991 

-0.7981 
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